Magnetic ordering in a doped frustrated spin-Peierls system 
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Based on a model of a quasi-one dimensional spin-Peierls system doped with non-magnetic im- 
purities, an effective two-dimensional Hamiltonian of randomly distributed S=l/2 spins interacting 
via long-range pair-wise interaction is studied using a stochastic series expansion quantum Monte 
Carlo method. The susceptibility shows Curie-like behavior at the lowest temperatures reached 
although the staggered magnetisation is found to be finite for T — > 0. The doping dependance of 
the corresponding three-dimensional Neel temperature is also computed. 

PACS numbers: 75.10.-b 71.27. +a 75.50.Ee 75.40.Mg 
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Quasi one-dimensional (ID) quantum antiferromag- 
nets exhibit fascinating magnetic properties at low tem- 
peratures. Some inorganic compounds, such as the ger- 
manate oxide CuGeOa* and the vanadate oxide L1V2O5 2 
are excellent realizations of weakly interacting frustrated 
spin- 1/2 chains. A spin-Peierls (SP) transition to a 
gapped dimerised ground state (GS) has been seen ex- 
perimentally in CuGeOsi, and theoretical calculations 3 
point toward a similar scenario in L1V2O5. Doping 
with non-magnetic dopants is realized experimentally in 
CuGeOa by substituting a small fraction of copper atoms 
by zinc (or magnesium) atoms*. An intriguing low- 
temperature phase where antiferromagnetism coexists 
with the SP dimerisation was observed*. Similar doping- 
induced antifcrromagnetic (AF) ordering has also been 
observed in the interacting dimer compound TlCuCbj 5 . 
Theoretically, each dopant is expected to release a soliton 
which can be viewed as a single unpaired spin separat- 
ing two dimer configurations*, hence leading to a rapid 
suppression of the spin gap under doping-. In an ideal- 
ized spontaneously dimerised spin chain the soliton would 
not be bound to the dopant-. The physical picture is in 
fact completely different: A static bulk dimerisation is en- 
forced by, e.g., couplings to the three-dimensional (3D) 
lattice, thus generating an attractive potential between 
the soliton and the dopan1p*-. 

The dopant-soliton confinement mechanism is respon- 
sible for the formation of local S=l/2 magnetic mo- 
ments*!*. Within a realistic model including an elas- 
tic coupling to a two-dimensional (2D) lattice*!, it 
was shown that these effective spins experience a non- 
frustrated interaction that could lead at T = to a finite 
staggered magnetization**. Recently, similar conclusions 
were reached using a model with purely magnetic inter- 
actions including a four-spin exchange coupling 12 . In this 
Letter, we analyze the formation of the dopant-induced 
AF order which coexists with the SP dimerisation. After 
using exact diagonalisation (ED) of small clusters (fol- 
lowing Refs»*i*i*) to construct an effective diluted S=l/2 
model, we take advantage of the non-frustrated char- 
acter of the resulting Hamiltonian to perform extensive 



state-of-the-art stochastic series expansion (SSE) quan- 
tum Monte Carlo (QMC) simulations on 2D lattices as 
large as 288 x 288 with up to N s = 576 (dopant) spins and 
down to temperature as low as T = 1/ j3 = 2~ 14 (or 2~ 18 
for N s — 256). The uniform susceptibility is shown to 
exhibit a Curie-like behavior although the T = and in- 
finite size extrapolated staggered magnetization is found 
to be finite down to the smallest dopant concentrations 
x available. The Neel temperature (assuming a small 
3D coupling) is also computed versus x and compared to 
experiments. 

We start with the microscopic Hamiltonian of a 2D 
array of coupled frustrated spin-i chains and we sum- 
marize the procedure followed in refs^*** to derive an 
effective Hamiltonian. 

H = ^^[7(1 + 6i >a ) Si <a ■ Si+l l0 

i,a 

+ CtJ Si, a ■ S J+2 , a + hi ia Sf ta ] , (1) 

where the i and a indices label the L sites and M chains 
respectively. The energy scale is set by the exchange cou- 
pling along the chain (J = 1) and a is the relative mag- 
nitude of the next nearest neighbor frustrating magnetic 
coupling. Dopants are simply described as randomly lo- 
cated inert sites (i, a) (see Fig. where S^ a = is set 
in Eq. (T). Small inter-chain couplings are included in a 
mean-field self-consistent treatment assuming, 



hi,a = J±{(Si. a +l) + ( S la-l)) 



(2) 



5i,a — ■j - {\Si )0 +i • Sj_|_i )a _|_i)-|-(S,- i0 _i • S,-_(_i )( j_i)}(3) 

While the first term accounts for first-order effects in 
the inter-chain magnetic coupling J±, the second term 
might have multiple origins; although a four-spin cyclic 
exchange mechanism provides the most straightforward 
derivation of i1i*£, at a qualitative level, J4 can also mimic 
higher-order effects in J_l 13 or the coupling to a 2D (or 
3D) lattice. In that case, due to a magneto-elastic cou- 
pling, the modulations 5i >a result from small displace- 
ments of the ions. The elastic energy is the sum of a 




local term \K\\ a 5f a and an inter-chain contribution 
^J- J2i a h,ah,a+i of electrostatic origin— In that case 
Eq. Q is replaced byJ£, K\\6 ita + K±(S i>a+1 + 5 i)0 _i) = 
J(Sj. a • Si+i )0 ), giving very similar results^ so that we 
shall restrict to Eq. © here. 

By breaking a dimer, each dopant releases a soliton 
carrying a spin 1/2. As shown previousl}^, J4 (or equiv- 
alently K± in the alternative model 1 -) leads to a con- 
finement of the moment to the dopant with a localiza- 
tion length ~ J^ v ■ At temperatures lower than the 
spin gap, these effective spins dominate the physics (see 
Fig^l and a low-temperature description of the doped 
system can be obtained using an effective model includ- 
ing only (long-range) pair-wise interaction between these 
local moments. The coupling J between two effective 
spins at arbitrary relative distance is computed by ED 
as the energy difference between their singlet and triplet 
configurationsii*i£. The sign of this interaction (i.e. its 
ferromagnetic or AF nature) depends on whether the two 
dopants lie on the same or opposite sublattices^, so that 
an overall AF ordering is favored. 

In order to derive an analytic expression for J cff valid 
at long distances, we fit the numerical ED data (re- 
stricted to short and intermediate length scales). A long- 
range non-frustrated Heisenberg model of diluted effec- 
tive spin-i can then be defined, 

U cS = e ri e r2 J cff (ri - r 2 )S ri • S r2 , (4) 

ri,r 2 

with e r = 1 (0) with probability x (1 — x), where x is 
the dopant concentration. Such a model can be stud- 
ied by QMC on L x x L y clusters much larger than 
those accessible to ED— and at all temperatures. Us- 
ing five phenomenological parameters, two energy scales 
and three length scales, simple expressions fit the ED 
data for a wide range of (physical) parameters. Nat- 
urally, one has to distinguish four cases depending on 
whether the dopants arc located on the same (Aa = 0) 
or on different chains, and whether they are located on 
the same or on different sublattices. When Aa = 
(same chain), J cff approximately fulfills </ off (At,0) = 



FIG. 2: Staggered magnetic structure factor per site vs 
inverse-temperature (3 computed using a /3-doubling scheme 
and averaged over 1000 to 2000 samples. Results for N s = 256 
spins on 56 x 56, 64 x 64, 72 x 72, 80 x 80, 88 x 88 and 96 x 96 
(from top to bottom) lattices correspond to the concentra- 
tions x indicated on the plot. 

Jo(l — Ai/^jj 1 ) for At even and At < £j? and other- 
wise J (At, 0) = 0. For dopants located on different 
chains and on the same sublattice (At + Aa even) one 
has, J oS (Ai,Aa) = - J' Q exp(-|i) exp(-|^) , while, if 

the dopants are on opposite sublattices, J cff (At,Aa) = 
J^exp(-As) for At < 2£|| and J off (Ai,Aa) = 

4 exp(- A '~ 2e " ) exp(-A2:) ; for At > 2^. The fitting pa- 
rameters used here for a = 0.5, J± = 0.1 and J4 = 0.08 
are J = 0.52, J' a = 0.3, ^ = 17.3, £|| = 2.5 and £_l = 1. 
We stress that the alternating sign of J is a crucial fea- 
ture of the interaction which guarantees the absence of 
frustration. Note also that the magnitude of J cS shows 
a unique exponential behavior oc exp(— A^ — As) at long 

distance although its short distance behavior is more 
complicated (but probably not relevant). The distribu- 
tion of these coupling contains a very large density of 
couplings of small magnitudes. 

We study the effective Heisenberg model using the SSE 
methodic to investigate GS as well as finite T properties. 
In this approach, the interactions are sampled stochasti- 
cally, and for a long-ranged interaction the computational 
effort is then reduced from ~ N% to N s In (jV s )— In or- 
der to accelerate the convergence of the simulations at the 
very low temperatures needed to study the ground state, 
we use a /^-doubling scheme^ where the inverse temper- 
ature is successively increased by a factor 2. Comparing 
results at several j3 — 2", one can subsequently check 
that the T — > limit has been reached. 

The AF ordering instability is signalled by the diver- 
gence with system size of the staggered structure factor, 

S{^) = -^—{C£{-l) i Stf). (5) 
Note that within our effective model approach, only the 




FIG. 3: Staggered magnetization per site. Disorder average 
has been done over at least 2000 samples, (a) Finite size 
extrapolations (see text) at fixed doping x. The symbols are 
identical to Fig|5]for x varying from ~ 2.8% to ~ 8.2%. (b) 
Doping dependance of ttiaf for various numbers of spins and 
in the thermodynamic limit (full symbols). 



sites carrying a "dopant spin" contribute to this sum. It 
is convenient to normalized S with respect to the num- 
ber of sites, i.e. to define a staggered structure factor 
per site; s(tt,tt) = S(tt, ir)/L x L y . In an ordered AF 
state, s(ir, n) should converge, with increasing size, to 
a non-zero value < 1/4. The (finite size) sublattice 
magnetization toaf can then be obtained by averaging 
s(tt, 7t) over a large number of dopant distributions, i.e. 
(ttiaf) 2 = 3(s(7r, 7r))dis, where the factor 3 comes from 
the spin-rotational invariance22> and (. . .)di s stands for 
the disorder average. The staggered magnetization per 
dopant is then simply TO sp i n = tuaf/x. We have checked 
that extrapolations to the thermodynamic limit using dif- 
ferent aspect ratios L y /L x give similar results and, here, 
we only report data for L y = L x = L. 

Since, strictly speaking, in 2D the divergence of S(tt, tt) 
occurs only at T = (f3 = oo) it is appropriate to first 
extrapolate the finite size numerical data to T = 0. As 
shown in Fig. |2] the staggered structure factor saturates 
at sufficiently low temperature and the GS value of ttiaf 
(averaged over disorder) can be safely obtained. Then, 
using a polynomial fit in \/\fW s (order 2 is sufficient) 
an accurate extrapolation to the thermodynamic limit, 
N s — ► co (or L — » oo at constant x), is performed as 
shown in Fig|3fa). The doping dependance of the extrap- 
olated rriAF is given in Fig|3fb). Note that our results, 
although consistent with previous (T = 0) extrapolations 
attempted on small clustersAi, are far more accurate due 
to the use of much larger systems. We have tested various 
fits to the data. Assuming a power law oc x^ , the best fit 
(solid line in FigE^b)) gives an exponent fi ~ 1.38 > 1. 
However, the alternative behavior a\x + a2X 2 would only 
be distinguishable at even smaller x. 

We have computed the uniform susceptibility x(T) for 



a wide range of temperatures. Results for the Curie con- 
stant C(T) = T\{T) are shown in Fig.H At the highest 
temperatures the effective dopant spins behave as free 
spins while at low temperature we observe a new Curie- 
like behavior with a reduced Curie constant ~ 1/12. Al- 
though here the 2D system orders at T = (as proven 
above) this behavior agrees with a qualitative argument 
by Sigrist and Furusak: 16 based on the formation of large 
spin clusters. A detailled analysis of a low-temperature 
scaling regime similar to the one observed in random 
ferromagnetic-antiferromagnetic spin chains^ will be re- 
ported elsewhere^. 

We finish this investigation by calculating the Neel 
temperature, asssuming a small (effective) 3D magnetic 
coupling A3D between the 2D planes. Using an RPA cri- 
terion, the critical temperature Tn is simply given by 
Xstag(? 1 A r ) = 1/|A3d| where the staggered spin suscepti- 
bility (normalized per site) is defined as usual by, 

X et jT) = f driSmsm , (6) 

i,j 

and averaged over several disorder configurations (typi- 
cally 2000). Since Xstag(Tjv) is expected to reach its ther- 
modynamic limit for a finite linear size L as long as T)v 
remains finite, accurate values of Tjv can be obtained us- 
ing a finite size computation of Xstag(^) for not too small 
inter-chain couplings. Fig. 0a) shows that Xstag(r) di- 
verges when T ~ * 0. Tn is determined by the intersection 
of the curve XstagfT) with an horizontal line at coordinate 
I/A3D. Note that finite size corrections remain small, 
even in the worst case corresponding to very small X30 
values and large dopant concentrations. The doping de- 
pendance of Tn is plotted in Fig. |5|b) for a particularly 
small value A3D = 0.01 (in order to show the small size 
dependance observable in that case). It clearly reveals a 
rapid decrease of Tn when x — > 0, but, in agreement with 
experiments, does not suggest a non-zero critical concen- 
tration. In FigEfb), wc show the behavior of Tn(x) down 
to x ~ 0.007. Note that from numerical fits of our data, 
we can not clearly distinguish between a power-law be- 
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havior (with an exponent ~ 2.5) and an exponential law 
like Aexp(—B/x), as suggested by fits of experimental 
data for Cui-^Zn^GeO^. 



In summary, we have studied an effective low-energy 
Hamiltonian (valid for low temperature physics) describ- 
ing the interaction of a finite concentration of spin- 
less dopants randomly distributed in a generic low- 
dimensional (frustrated) SP system. The SSE method, 
which is applicable here because of the non- frustrated na- 
ture of the effective model, was used in combination with 
finite-size scaling to compute both GS and finite temper- 
ature properties. The uniform susceptibility exhibits a 
Curie-like behavior down to very low temperature. We 
also predict that the AF order developps continuously 
without any finite critical dopant concentration in agree- 
ment with experiment. 
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FIG. 5: (a) Staggered susceptibility of a 2D layer vs temper- 
ature (using log-log scales) for N 3 = 256 (full symbols) and 
N s — 576 (open symbols) spins. Concentrations x are shown 
on the plot (similar symbols as in Figs. H and EJ. = 100 
is shown by the dashed line, (b) Neel temperature vs dopant 
concentration x for a 3D RPA inter-plane coupling X^d = 0.01 
and for N s — 256 and N s = 576 spins. 
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